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We show how to consistently construct initial conditions for the QCD evolution equations for 
double parton distribution functions in the pure gluon case. We use to momentum sum rule for 
this purpose and a specific form of the known single gluon distribution function in the MSTW 
parameterization. The resulting double gluon distribution satisfies exactly the momentum sum 
rule and is parameter free. We also study numerically its evolution with a hard scale and show 
the approximate factorization into product of two single gluon distributions at small values of x, 
whereas at large values of x the factorization is always violated in agreement with the sum rule. 
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I. INTRODUCTION 

Multiparton interactions play an important role in the hadronic collisions at high energies. They occur when at one 
encounter of the initial hadrons, more than one partonic interaction occurs. They were first observed and measured 
at the Tevatron 1-4] and subsequently a systematic experimental study was performed at the higher energy Large 
Hadron Collider @-u|. The theoretical description of such interactions within perturbative QCD is possible in the 
presence of the sufficiently hard scales. The computation of double parton scattering (DPS) cross sections within the 
collinear framework makes use of the double parton distribution functions (DPDFs) Ii 8j-j35l . In the collinear leading 
logarithmic approximation DPDFs obey QCD evolution equations [1, 9t, EH, EH, EH PIGI . similar to the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations for single parton distribution functions (PDFs). The 
evolution equations for DPDFs conserve new sum rules which relate the double and single parton distributions once 
they are imposed on initial conditions for the evolution equations at some initial scale. All the attempts up till now 
to construct conditions which satisfy these sum rules were rather unsuccessful, see e.g. Refs. EH, EH, EH] with an 
exception of the analysis [38] for valence quarks only. 

In this letter, we show how to consistently perform such a construction in a pure gluon case, using the known single 
PDFs in the MSTW parameterization [39] and the momentum sum rule. We find the parameter free double gluon 
distribution which we evolve with our numerical program. In particular, we study the build up of its approximately 
factorizable form for small values of parton momentum fractions, 2 : 1,2 < 0.1. The full case with quarks and gluons is 
postponed to a separate publication. 


II. EVOLUTION EQUATIONS AND SUM RULES 

We consider the DPDFs with equal hard scales, Qi = Q 2 = Q , and the relative momentum q = 0: 


D fl f 2 (x 1 ,x 2 ,Q) = D flh (x 1 ,x 2 ,Q,Q,q = 0) 


(1) 
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where £ 1.2 G [0,1] are parton momentum fractions, which obey the condition x± + X 2 < 1, and /i t 2 are parton flavors 
(including gluon) 17|,|22|. In this case, the evolution equations in the leading logarithmic approximation read 


din Q- 


D flh (x 1 ,x 2 ,Q) = 


du . 
u 


[ D flf ,{xi,u,Q) 

Jx 2 U \ U / 

—j— Pp^f, f 2 ( —X—) D fl (x 1 +x 2 ,Q)\. 
X1+X2 1 1 J \X 1 +x 2 J j 


( 2 ) 


where the functions V on the r.h.s. are the leading order Altarelli-Parisi splitting functions (with virtual corrections 
for Vff included). The third term on the r.h.s corresponds to the splitting of one parton into two daughter partons, 
described by the Altarelli-Parisi splitting function for real emission, Pp^f, f 2 ■ It contains the single PDFs, Df, thus 
eq. © has to be solved together with the ordinary DGLAP equations, see e.g. Ref. [Hj for more details. 

The significance of the splitting terms in the evolution equations © for the computation of the double parton 
scattering cross sections was a subject of intensive debate in the literature over the last few years [l7|, [3 0, HU, 0, 
[3ll . fibi . Hcj . The conclusion which emerges from this discussion is that the processes which are summed up by the 
splitting terms and coming from both hadrons in hadron-hadron collisions should rather be classified as the single 
parton scattering process [26;]. On the other hand, the so called single splitting contributions, with parton splitting 
from one hadronic side only, are important for the double parton scattering cross sections [3, HI HI, El ■ From the 
perspective of the present paper, in which we only concentrate on the evolution of the DPDFs, the splitting terms in 
the evolution equations are crucial for the conservation of sum rules which are discussed below. 

The sum rules which are conserved by the evolution equations © are the momentum and valence quark number 
sum rules El- Imposing them for initial conditions specified at some initial scale Qo, they are guaranteed to be 
satisfied at any other scale Q. The momentum sum rule for the DPDFs reads 

r 1 — 3,2 

E/ dx 1 x 1 D fl f 2 (x 1 ,x 2 ) = (l-x 2 )Df 2 (x 2 ) , (3) 

, Jo 


while the valence quark number sum rule is given by 


1 —X 2 

dx 1 {D qf2 (x 1 ,x 2 ) - Dqf 2 (xi,x 2 )} = (N t j - S hq + S hq )D f2 {x 2 ), (4) 

where q = u,d, s and N u = 2, Nd = 1, N s = 0 are the valence quark number for each of the quark flavors. The same 
relations hold true with respect to the second parton 



dx 2 X 2 D h f 2 {x l ,X 2 ) = (1 - Xi)D fl (x{) 


r*l— x\ 


dx 2 {D fiq ( xi,x 2 ) - Df iq (xi,x 2 )} 


{Nq Sf iq + 5f iq )D f x {x 1 ). 


Notice that if the DPDFs are parton exchange symmetric, 


( 5 ) 

( 6 ) 


D f 1 f 2 (xi,x 2 ) = D hfl {x 2 ,x 1 ), 


( 7 ) 


the sum rules with respect to the first parton imply the sum rules with respect to the second one since the evolution 
equations also conserve parton exchange symmetry. 

We see that the above sum rules relate the double and single parton distribution functions, which reflects the 
common origin of those distributions, namely, the expansion of the nucleon state in Fock light-cone components El- 
In addition, the sum rules for the single parton distributions are also satisfied - the momentum sum rule 


E 

/ 


/ dxxDf{x) = 1 
Jo 


( 8 ) 


and the quark valence sum rule for q = u,d : s 



{D q {x) - D q {x)} = N q . 


(9) 
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III. MELLIN MOMENT FORMULATION 


Let us perform the double Mellin transform of the DPDFs 

D f if 2 (ni,n 2 ) = f dx 1 [ dx 2 {xi) ni ~ 1 {x 2 ) n2 ~ l D flh (®i, x 2 )0(l - xi - x 2 ). (10) 


Jo Jo 

where ^ are complex numbers and we omit the scale Q g in the notation from now on. The step function 0(1— x\ —x 2 ) 
is inserted into the definition of the Mellin transform since this is the region over which the double parton distribution 
is defined. Similarly, for the single parton distribution functions, we define the Mellin moments 

r i 


Df(n)= / dxx n 1 Df(x ), (11) 

Jo 

where n is a complex number. The Mellin moments can be transformed back to the s-space using the inverse 
transformation for the single parton distribution, 


Df {xi) = J c ( x i) n Df{n ), 

and similarly for the double parton distribution function 

D flh (xi,x 2 ) = j c \ Oi ) - " 1 Jc 7^ (x 2 )~ n2 Df lh {n i,n 2 ), 


( 12 ) 


(13) 


where the integration contours C\ and C 2 lie to the right of the rightmost singularity in the complex plane of n\ and 
n 2 , respectively. Let us emphasize that formula m is only applicable to X\. 2 £ [0,1] and X\ + x 2 < 1 . 

The sum rules © and © can be written with the help of the Mellin moments after the integration of both sides 
over x 2 with the factor (cc 2 )" 2-1 . Thus, we find 


h 

= D h [n 2 ) - Df 2 (n 2 + 1), 

(14) 

D q f 2 {l,n 2 ) - Dq h (l,n 2 ) 

= ( N q — Sf 2 q + 5f 2 q)Df 2 (n 2 ). 

(15) 

Analogous relations hold true for the second parton 



^ Afi/ 2 ( n i> 2 ) 
h 

= D h {m) - D h {n 1 + 1), 

(16) 

Dfiq{ni,l) ~D fl q(ni,l) 

= {.Nq — Sf iq + 5f 1 q)Df 1 {n\). 

(17) 


These sum rules have to be satisfied simultaneously with the momentum sum rule for the single parton distribution 


J2D f (2) = l, (18) 


and the valence quark sum rule 


D q { 1) - Af( 1) = N q . (19) 

It would be extremely useful to construct initial conditions for DPDFs which fulfill the above sum rules since 
the PDFs on the r.h.s of Eqs. ©-© are very well known from the global analysis fits. Thus, the PDFs constrain 
the DPDFs, solving or significantly reducing the problem of uncertainty in the specification of initial conditions for 
DPDFs evolution. For this purpose, we consider the single PDF parametrization from the MSTW fits |39|. We will 
choose the leading order (LO) version of this parameterization since the evolution equations © are given in leading 
logarithmic approximation. We start from considering the case with only gluons. This limits the set of the possible 
distributions in which D gg and D g are the only relevant functions and of course we only have to fulfill the momentum 
sum rule. 
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IV. PURE GLUON CASE 

The single gluon distribution is specified in the LO MSTW parameterization at the scale Q o = 1 GeV and is given 
in the form 


D g (x) = A g x 5g 1 (l-x)‘ na (l + e g y/x + % g x), (20) 

where A g = 0.0012216, S g = —0.83657, = 2.3882, e g = —38.997 and = 1445.5. For errors on these parameters 

and the discussion of their determination, see (393- Since we only use gluons in our analysis, we renormalize the gluon 
distribution such that the total longitudinal momentum carried by gluons equals one, which results in A g = 0.0033723. 
This is really not so crucial here as the normalization can be set arbitrarily for the single channel case and it does not 
affect the subsequent discussion. The parametrization (1201) can be written in a general form which is more suitable 
for our purpose 

L 

D g {a ;) = XX**“2(1-*)#, (21) 

fc=l 

where L = 3 and the parameters N£, a k and /3 g can easily be found by the comparison with eq. (BUI) . 

N g = A g , Ng = e g A g , N 3 = 7 g A g 

a 3 = s 9 - 1 > a l = s 9 - 3 . a 3 g = 6 g , = P 2 g = (3 3 = rig. ( 22 ) 

In the Mellin space, the gluon distribution (1201) can be written as 


D g {n) 


r (n + a g fc )T(^ + l) 
“ 9 r(n + a k g+/3$ + 1) ' 


(23) 


where the expresion on the r.h.s., T(a;)r(?/)/r(a; + y) = B(x,y), is the Euler Beta function. Thus the MSTW 
parametrization for the initial condition is in the form of the sum over the Beta functions with different sets of 
parameters which govern the small x —> 0 and large x —> 1 behavior. 

For the double parton distribution D gg we shall take the following ansatz in the form 

L 

Dgg(x i,x 2 ) = ^Ngg {xix 2 ) aa (1 - Xi - X2Y 9 , (24) 

k =1 


where N k g ,a k and /3 g are the parameters to be determined. The above ansatz is in the form of the sum over the 
Dirichlet-type distributions of order K = 3 


if =3 

/(a:i,£2;71,72,73) = N n X T . ( 25 ) 

i=l 

with X\,X 2 > 0,xi + X 2 < 1 and 23 = 1 — x\ — 22 . Notice that the function (1241) is parton exchange symmetric, 
D gg { 21 , 22 ) = D gg (x 2 ,xi), as it should be. It should also fulfill the momentum sum rules with respect to both partons. 
Also note that the number of terms in this sum, L , is the same as the number of terms in the single parton distribution 
m ■ The Mellin space representation of the above ansatz reads 


L 


Dgg{n 1 ,n 2 ) = J2 N gg 
k= 1 


FK + dg)r(n 2 + Qg)r(l + p k g ) 

r(m + n 2 + 1 + fig + 2a k ) 


(26) 


which is in the form of the generalized Beta function. 

In the pure gluon case only the momentum sum rule for the DPDFs in the momentum space reads, 

1 — X 2 

dx 1 x 1 D gg (x 1 ,x 2 ) = (1 - x 2 )Dg(x 2 ), (27) 

and similarly for the momentum sum rule with respect to the second gluon. In the Mellin representation this condition 
reduces to 


Dgg(2,n 2 ) = D g (n 2 ) - D g (n 2 + 1 ), 


( 28 ) 
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It is easy to see that the distributions of the form presented in Eqs. m and (1241) fulfill the momentum sum rule 
provided certain constraints are satisfied. The right hand side of Eq. (1281) is the difference of the moments of single 
parton distributions which can be written as 


D g (n 2 ) - D g (n 2 + 1) = 


L 

E 

k=1 


N k 


B(n 2 


,#+ 2 )= 5> 4 


k =1 


k T(n 2 + a k g )T(l3 k g +2) 
9 T (n 2 + a k g +M + 2) 


(29) 


where we used the following property of the Beta function 


B{a,b) = B(a+ 1,6) + B(a,b+ 1) . (30) 

On the other hand the left-hand-side of Eq. (051) is obtained upon setting m = 2 in Eq. (HO) 


D gg %n*)=Y R i 


k 

99 ~ 


T(2 + a k )T(n 2 +d k )T(l + p k ) 


k=1 


T(3 + n 2 + P k q +2a k ) 


(31) 


Now in order for the momentum sum rule to be satisfied, we need Eqs. (1291) and (1311) to be equal term by term in the 
sum over k. From the requirement that the poles and zeros in n in each term should be in the same location we find 
that 


a k =a k , 2a k g +p k g +3 = a k +/3 k g +2, (32) 

and from the requirement that the normalization of each terms should be the same we have that 

N k g r(2 + a k )T(l + p k ) = N k T(2 + (3 k ). (33) 

From these relations we compute all the parameters of the double gluon distribution in terms of the known parameters 
of the single gluon distribution, given by Eq. (l22l) , to find the following parameter-free double distribution at the initial 
scale Qo = 1 GeV, 


Dgg(X 1 ,X 2 ) = ^ 


N„ 


nft + 2 ) 


k =1 


T(a k + 2)T(p k - a k ) 


{xix 2 ) a n (1 - Xi - a: 2 ) /3s “s’ 


(34) 


satisfing the momentum sum rule (1271) by construction. Notice that even for small momentum fractions, x\ t2 1, 
the resulting double gluon distribution is not factorizable, i.e. D gg (x\,x 2 ) 7 ^ D g {x\)D g {x 2 ). 


V. EVOLUTION OF DOUBLE GLUON DISTRIBUTION 


The evolution equations ([2]) reduced to the pure gluon case have the following form 

rl2 du 


d ri < as 

—2 D gg (x 1 ,x 2 ,Q) = 


din Q 


^ ^ B gg (u,x 2 , Q) 

f x ~ Xl du (X2 \ . 

+ / —V gg (—JD gg (xi,u,Q) 

Jxn u \ u / 


+ 


1 


X\ + x 2 


■ p. 


gg 


Xl 


X\ + x 2 


D g (xi +x 2 ,Q) >. 


(35) 


where is the gluon-to-two gluon splitting function for real emission in the LO approximation. Strictly speaking, 
such an equation can be reasonable approximation for small values of momentum fractions. Using our numerical 
program, we solve the above equation with the initial condition (1341) . We compare our results with those obtained 
from the usually assumed form of the initial conditions [l 8 | , which satisfy the momentum sum rule only approximately, 


where the correlation factor 


Dgg(x l,x 2 ) = Dg(Xl)Dg(x 2 )p(Xl,X 2 ) 
(1 - Xi - x 2 ) 2 


p(xi,x 2 ) = 


(1 - a:i) 2 (l - x 2 ) 2 


(36) 


(37) 
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x 2 =0.01 




x i 




FIG. 1: The distribution xiX 2 D gg (xi, xi = 10 2 ) at Qq = 1 GeV 2 (left upper panel) and Q 2 = 10 GeV 2 (right upper panel) 
and the ratio (13811 (lower panels). The solid lines correspond to input (1341) (our) while the dashed lines to input (1361) (prod). 


The results are shown in Fig. |T] and Fig. [2] We plot there the double gluon distribution X\X 2 D gg (x\, X 2 ) as a 
function of x\ for two values of the scale, initial Qq = 1 GeV 2 and Q 2 = 10 GeV 2 (upper panels), for two fixed fixed 
values of x 2 , small 10 -2 and large 0.5, respectively, The solid lines show the results obtained from our input (l34l) while 
the dashed lines correspond to the input (l36l) with the gluon distribution (l20l) . In the lower panels we plot the ratio, 


ratio = 


Dgg(x u x 2 ) 
Dg{xi)Dg{x 2 ) ’ 


(38) 


which characterizes factorizability of the double gluon distribution into a product of two single gluon distributions. 

For both values of x 2 , the initial double gluon distributions differ significantly for small values of x\, up to 10 -1 for 
X 2 = 1CV 2 and up to 10~ 3 for X 2 = 0.5. However, the QCD evolution equation erases this difference already at the 
scale Q 2 = 10 GeV 2 , see the upper panels in both figures. As we have observed, the initial distribution (ITU) is not 
factorizable into a product of two single gluon distributions for any values of x\ and X 2 - However, if both momentum 
fractions are small (< 0.01,) D gg becomes factorizable with good accuracy after evolution to the shown value of Q 2 , 
see the lower panels in both figures. A small breaking of the factorization can be attributed to the non-homogeneous 
term in the evolution equation (1351) . If one of the two momentum fractions is large, like the shown X 2 = 0.5, this is 
no longer the case and the factorization is significantly broken for all values of x± independent of the values of the 
evolution scale. We check that for larger values of Q 2 that than shown here. We have to remember, however, that 
the large x domain has to be supplemented by quarks. 

In conclusion, the initial double gluon distribution (1341) is very different from that proposed so far. However, the 
QCD evolution equation significantly diminishes this difference at not so high values of the evolution scale Q 2 . 


VI. SUMMARY 

In this paper we constructed the double gluon distribution D gg from the known single gluon distribution D g , given 
by the MSTW parameterization, in the pure gluon case. The construction is based on the expansion in terms of 
functions which are the Dirichlet distributions. Since the MSTW distribution has already the form of the sum over 
Dirichlet distributions of order 2, we postulated the double gluon distribution as a sum over the Dirichlet distributions 
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FIG. 2: The distribution xiX 2 D gg (xi, £2 = 0.5) at Q 2 = 1 GeV 2 (left upper panel) and Q 2 = 10 GeV 2 (right upper panel) and 
the ratio (13811 (lower panels). The solid lines correspond to input (1341) (our) while the dashed lines to input (1361) (prod). 


of order 3 and identified the parameters in the two distributions using the momentum sum rule for this purpose. As 
a result, we obtained the parameter free double gluon distribution which we evolve with the QCD evolution equation. 
We studied the build up of the factorized form of D gg with the increasing evolution scale Q. We found that such a 
form approximately sets up rather quickly for small momentum fractions, 2 : 1,2 < 0.1. As expected, for higher values 
of x, the factorized form is not valid at all due to the constraint from the momentum sum rule. 

The next step would be to extend this formalism to include the quarks and satisfy the momentum and valence 
quark sum rules simultaneously. The expansion in terms of the Diriclilet functions can be constructed also in the case 
with quarks. The whole formalism is however much more complicated due to the large number of the double parton 
distribution and therefore the full analysis will be presented in the future publication. 
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